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Faddeev and Jackiw’s method for constrained systems is used to derive a gauge invariant formula¬ 
tion of cosmological perturbations in the one bubble inflationary universe. For scalar perturbations 
in a flat universe, reduction of the action to the one with a single physical degree of freedom has 
been derived in the literature. A straightforward generalization of it to the case of an open universe 
is possible but it is not adequate for quantizing perturbations in the one bubble universe, because 
of the lack of Cauchy surfaces inside the bubble. Therefore we perform the reduction of the ac¬ 
tion outside the lightcone emanating from the center of the bubble or nucleation event, where the 
natural time constant hypersurfaces are no longer homogeneous and isotropic and as a result the 
conventional classihcation of perturbations in terms of scalar and tensor modes is not possible. Nev¬ 
ertheless, after reduction of the action we find three decoupled actions for three independent degrees 
of freedom, one of which corresponds to the scalar mode and the other two to the tensor modes. 
Implications for the one bubble open inflationary models are briefly discussed. As an application 
of our formalism, the spectrum of long wavelength gravity waves is simply obtained in terms of the 
real part of the reflection amplitude for a one dimensional scattering problem, where the potential 
barrier is given in terms of the bubble prohle. 


I. INTRODUCTION 

Models which reconcile inflation with a non-critical density, O 1, have been recently proposed in the literature 
[Q. Although these models are somewhat more involved than standard inflation, they may end up being favored by 
observations ||J^. In this scenario, one starts with a scalar held trapped in a false vacuum that drives a de Sitter-like 
phase of inflation. This held undergoes a first order phase transition by forming an 0(3,1) symmetric bubble. Inside 
the bubble, the scalar held in the new phase slowly rolls down the potential, driving a short second period of inflation. 
Our observable universe would be contained inside a single bubble, whose symmetry accounts for the observed large 
scale homogeneity and isotropy. The second period of inflation is needed to solve the entropy problem. 

A complete study of cosmological perturbations in open inflation involves the quantization of fields in the presence 
of a bubble. So far, progress has been made by quantizing the scalar field but ignoring the selfgravity of these 
fluctuations. 0-^. The quantization of tensor modes has been considered in Some interesting features have 

been found. There are some scalar modes - the so called supercurvature modes- which are not normalizable in the 
open hyperboloids but do contribute to the microwave background anisotropies. Some of these, which correspond to 
fluctuations of the bubble wall, are found to be such that they can be rewritten as tensor modes, due to their special 
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eigenvalue of the Laplacian 110. Of course, it would be very interesting to know whether these features survive 
when gravitational perturbations are included and all the perturbations are consider in a unified manner. 

Gravitational perturbations contain gauge degrees of freedom. Perhaps the most elegant way to get rid of the 
“unphysical” gauge modes is the gauge invariant theory of cosmological perturbations. |Iq-El|. In a universe domi¬ 
nated by a single scalar field, there is only one physical degree of freedom for scalar perturbations. Hence to quantize 
perturbations, one has to reduce the number of variables in the action to a single variable by using the constraints. 
This program has been carried out for scalar perturbations in a spatially flat universe and in a spatially closed 

universe |24|] . The final form of the action resembles the one for a scalar field with a time-dependent mass term. One 
might expect that the extension of this program to the case of an open universe would be straightforward. In fact, it 
would be so if the hypersurfaces of homogeneity and isotropy of an open universe were Cauchy surfaces on which the 
canonical commutation relations could be set up. However, in the case of one bubble inflationary scenario, the whole 
universe is contained in a single bubble and the open hypersurfaces foliate only the interior of the lightcone emanating 
from the center of the bubble [0. In particular, one cannot deal with supercurvature modes on these hypersurfaces. 
Hence the quantization should oe carried out outside the lightcone where Cauchy surfaces exist but no hypersurface 
of homogeneity and isotropy exists. The purpose of this paper is to carry out this non-trivial task, i.e., to find the 
reduced action for cosmological perturbations appropriate for the scenario of the one bubble inflationary universe. 

Dirac’s procedure has been for a long time the canonical way to treat constrained systems. Faddeev and Jackiw 
(FJ), however, have proposed 0| an alternative approach which leads to the same results without following all of 
Dirac’s steps. As they point out, two aspects of Dirac’s procedure can be avoided. First, it is not necessary to 
distinguish between different classes of constraints: all of them can be treated on equal footing without ambiguities. 
Second, it is not necessary to define conjugate momenta for those velocities which appear linearly in the Lagrangian, 
as is customary done in Dirac’s approach. Applied to our case, the method gives the linearized action for cosmological 
perturbations in terms of three gauge invariant degrees of freedom, one corresponding to the scalar mode and two to 
the tensor modes. This action is then ready for canonical quantization. 

The paper is organized as follows. In Section H we describe the method of reduction. In Section HI we apply it to 
cosmological perturbations in open inflation. In Section IV we derive, as an aplication of our formalism, the spectrum 
of gravity waves in open inflation. In Section V we summarize our conclusions. Some technical issues are left to the 
appendices. 


II. REDUCTION METHOD 

In the Faddeev-Jackiw approach one begins with an action first order in time derivatives, 

S = J dt, (2.1) 

where are the phase space variables of the system, and z are a subset of those which do not appear in the kinetic 
term. The basis of the method 0] is to use the Euler-Lagrange equations of motion that contain no time derivatives 
(the real constraints of the theory) to reduce the phase space. The equations of motion of the z coordinates belong 
to this category. One starts by solving this set of equations first for as many z’s as possible and then, if there are 
any z’s that appear linearly in the Lagrangian, for as many ^’s as possible. After substituting these relations into the 
original action, it takes the form 


s* = I + L*{U)) dt, ( 2 . 2 ) 

where now the label i spans fewer coordinates than fj,. From now on, a * means that the known constraints have been 
substituted. 

If there are further constraints in the theory, they manifest themselves as c ombinations of the equations of motion 
which contain no time derivatives. Writing the equations of motion from ( |2.2| ) in the form 

dS* BT* 

= he + ^ - dth =■■ he + Gee, t) = 0, (2.3) 

where := difj — djfi, each zero mode of the kinetic matrix /y, i.e. ^Ih ~ us the constraint equation 

eiGeet) = o. ( 2 . 4 ) 

These constraints can be used again to reduce the phase space. The process is repeated until we end up with a nonsin¬ 
gular fij, which indicates that we have identified the reduced phase space in which the Lagrangian is unconstrained. 
The equivalence of this method with Dirac’s is discussed e.g. in P8| . 

For our present purpose, we do not have to follow the FaddeewHackiw approach step by step. In the problem of 
cosmological perturbations, the constraints are first class and become the generators of gauge transformations. In 


2 



this specific case, we can make use of this fact for reduction. As is well known, the first order quadratic action for 
gravitational perturbation takes the form. 


5 = J Par [p^, q^Q]) 


(2.5) 


where Cf^ are the gauge transformation generators, i.e. the gauge transformation of a variable is given by the Poisson 
bracket as Sgf = I f, J d^x >, and Q" = {Q’^^Pa) stands for the rest of canonical coordinates. The fields SAf'^ are 


the Lagrange multipliers corresponding to the perturbation of the lapse function and the shift vector. In our present 
problem, constraints Cg = 0, which are linear in the perturbed variables, are easily solved for as p^ = Pfj,[r jQ^]- 
That is, by taking a linear combination of C^, we can rewrite the constraints as — Pfj.[q'^,Q^]- Accordingly, 

we write Sgf = I f,J A'^ (7^ d^x >. Substituting p^ = p^[q‘', Q"] back into the action, we obtain 


( 2 . 6 ) 


s* = J {Mq^,Q^]r+Par-'H*)d^x, = 


which loses its canonical form but still keeps gauge invariance. This can be deduced from the well known fact that 
the constraints are first class: {0^,(5^ = 0.|^ Then we find 

dgPt. = {prr,Ql,J A“a d^x^ . (2.7) 

Now, we can show the gauge invariance of S* as 


SgS* = \p^[q'',Q% / A“C„d3a; 


SS 

Spfj. 


. a 


+ SgQ^ 


SS 


c„=o 


SQ^' 


Cu=0 


= SgS 


= 0 . 


( 2 . 8 ) 


c„=o 


Writing the gauge transformation of a given variable in terms of operators acting on the gauge parameters, i.e. 
SgQ'^ =■ SQ^[X^], gauge invariance can also be written as 


which implies that 






SS* ^i,SS*, 

sy = isyi' 


5* = 0, 


(2.9) 


( 2 . 10 ) 


where dQ" is the operator conjugate to SQ"^ defined as 


d'^xSQri [f]g:= / d'^x/SQ^g], 


( 2 . 11 ) 


for any square-integrable functions / and g. Thus, recalling the chain rule, S* should depend on q^ only through the 
following combination of variables 


Q" := Q" - =Qn_>QUl . 


( 2 . 12 ) 


**In general, this equality holds only in a weak sense. But, in our present problem, Cg are linear in the variables and hence 
their Poisson brackets are just numbers. Therefore this weak equality can be interpreted as a strong one. 
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It is easily shown that the new set of reduced variables Q” is gauge invariant. In an actual calculation, all we have to 
do is just to set = 0 in S* and reinterpret Q" as their gauge invariant counterparts, Q". In this way, we go from 
the set of va riable s {Q^,q^,P)j.) to the reduced set Q". 

Equation (2.10|) can be also derived directly from the full Faddeev-Jackiw procedure. If instead of taking advantage 


of the gauge invariance we followed the method step by step, after subst itut ing the four constraints we would find 
that the kinetic matrix has four zero modes. The constraint equations for these zero modes turn out to be 


SS* 

^+^Qu 


.5S* . 


= 0 , 


(2.13) 


but, as we have seen, due to the gauge invariance, the left-hand side vanishes identically. As explai ned, t hese identities 
point out that the action can be written in terms of the reduced gauge invariant set of variables (2.12). 


III. OPEN INFLATION 


As mentioned in the introduction, the interior of an 0(3,1) symmetric bubble is isometric to an open Friedmann- 
Robertson-Walker (FRW) universe. Unfortunately, the hypersurfaces of homogeneity and isotropy of this universe are 
not appropriate for setting canonical commutation relations or normalizing modes [ P5| , because they are not Cauchy 
surfaces for the whole space time. Therefore, we shall need to quantize on spacelike hypersurfaces which cut right 
through the bubble, and which are therefore not homogeneous. This renders the decomposition into scalar, vector and 
tensor modes into a somewhat unfamiliar form. In the end, however, the three standard physical degrees of freedom 
will be identified. 

The open FRW chart 

ds^ = —dt^ + a{t)‘^dQH 3 , dflns = dr^ + sinh^ r{dd'^ + sin^ dp'^), (3.1) 

covers only the interior of the lightcone emanating from the center of the bubble t = 0, r = 0, which we shall call the 

nucleation event N. Following |^, we shall call the interior of the lightcone region R and the outside of this lightcone 
region C. Region C can be covered by analytically continuing the coordinates t and r in region R to the complex 

plane. By taking t = ir and r = x -I- i(7r/2), with r and x real, the line element becomes 

ds^ = dr^ + aE{T)‘^dflds, dflds = "/ffdx^dx^ = —dx^ + cosh^ xi^O^ + sin^ dip'^), (3.2) 

where a£:(r) = —ia(ir), and dQdS is the metric of a (2-1-1) dimensional de Sitter space. In this chart, t is a ‘radial’ 
spacelike coordinate, whereas x is timelike. Now the spacelike hypersurface x = 0 is a Cauchy surface for the entire 
space-time p^ . It is convenient to introduce the conformal ‘radial’ coordinate rjE, with dr = —aEdrjE- Close to the 
lightcone emanating from the nucleation event N, the scale factor behaves as aEij) « r, and r]E -l-oo. As we move 
away from N along the r direction, the scale factor rises to a maximum and then decreases again, reaching another 
zero at the so called antipodal point A, which corresponds to r]E —oo. 

Writing the perturbed line element and the perturbed scalar field in the form 

= aEi'nE)‘^{{'^ + 2A)dr]l, - 2Sidx'-dr]E + ixff - hij)dx'-dx^}, , . 

P = Po{ve) + d(/7. 


the second order action for small perturbations is given by 




I d^xa%y^{-2{2n^ + n')A^ + (%|,) - -ff - ( 5 ^ " 

+k{S(p''^ + - 2k{(PqS(p'A - a%V^6(pA) 

-4(5*p - - HA) - (2A2 + 2Ah\ - 


(3.4) 


where k = SttG, = a'/a and a prime denotes a derivative with respect to the conformal ‘radial’ coordinate tje- 
This can be found e.g. from Appendix ^ just replacing ~a|; and jij ~lif in expression (|A8|). 

As mentioned before, the usual expansion of the metric perturbations (see e.g. BI) in scalar, vector and tensor 
modes with respect to the 3-hyperboloid, used inside the lightcone, cannot be used outside. The reason is that, in 
region C, the corresponding 3-hyperboloid on which these harmonics of various types are defined no longer gives a 
spatial section of the spacetime. Instead we shall expand in scalar and vector modes with respect to the 2-sphere. 
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Using a conformal time-like coordinate p defined through the relation coshx dp = dx, the metric element for the de 
Sitter space can be written as 


dsis = csipfi-dp"^ + iUABdx^dx^), (3.5) 

with ce{p) = cscp, and where A, B, ... run over 9 and ip. For convenience, we define He = ceIce, where a dot 
indicates derivative with respect to the time-like variable p. 

We write the metric perturbations as 

Sa = T\\a + Va, 

hpp = 2c%C, (3.6) 

hpA = —(?e (C||^ + WOl), 

hAB = C% {{w + ^'^'^Av)uiab -^v^ab + “^FiAWB)), 


where ||yi and stand for the covariant derivative and scalar Laplacian associated with loab, respectively. The 

meaning of when operates on a vector or tensor quantity is explained immediately below. The fields S', T, C,, 
w and V are scalar modes, and Va, Wa and Fa are divergenceless vector modes with respect to the 2-sphere metric 
loab- More explicitly, we can write, say, S = J2im S^™(/9)F^™(U) for scalar modes and Va — (^) 

for vector modes, where y^'"(U) are the ordinary spherical harmonics which satisfy AF^"‘(U) = —l{l + l)F^’”(f2), 
EAB is the unit anti-symmetric tensor on the unit 2-sphere (egip = sin0 etc.) and = co^'^ecA- Now we can clearly 
state the meaning of It should be understood just as —l{l + 1) when it is decomposed into modes. In order to 
avoid writing the summation over I and m for notational simplicity, we use ^^^A instead of —I{I -I- 1). 

The scal ar and vector mo des transform differently under the parity transformation. Thus, inserting the decompo¬ 
sition (^) into the action (^), the scalar and vector modes decouple, so they evolve independently. According to 
the change of signature under the parity transformation, we refer to the scalar (vector) modes as even (odd) parity 
modes. 

We shall see that the odd parity modes contain one physical degree of freedom, which corresponds to odd parity 


tensor modes when analytically continued inside the lightcone [|10| , [ll[ . The even parity modes contain two degrees of 
freedom, one corresponding to the usual scalar and the other to even parity tensor modes. 


A. Odd parity modes 


First we consider the odd parity modes. The result provided in this subsection is essentially the same as that 
given in |pd]| . However, our present approach based on the FJ method is quite different from the conventional Dirac’s 
method used in the previous work pd| ]. Furthermore, the analysis of the even parity modes discussed in the next 
subsection is rather complicated compared with the odd parity modes. So, also to understand our strategy, it will be 
convenient to present the analysis of the odd parity modes first. 

The Lagrangian density for odd parity modes is 

[{Va - 2hEVA + clW'Af + (U^ - + 2)(U4 - cIF'a) 

-cl{F^ + IU^)((2)A + 2){Fa + Wa)] ■ (3.7) 


By the definition of conjugate momenta, we find 

:= ^^ = ^M^(v^-2hEV^ + clW^'), 

dVA 2 k.ce 

jjA a|ci;v/^((")A + 2) 

^ ■ dpA 2 k ^ 

Here we have raised indices with . If we cast it into first order form we obtain 


= n^U4 + n^FA, 

= -H^^^-- ^Ufa 


a%CE^/idjp^'^A + 2) 


kce 




n^nvA + 2hEii^VA 


(3.8) 


(3.9) 
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4 KCE 

Cit = -n^ + 4n^'. 


{V^ - + 2){Va - clF'^), 


To find the reduced phase space of the lagrangian (3.9), we solve = 0 for 11^, and substitute it back into the 
Lagrangian. Using t he fac t that Cw is the generator of odd parity gauge transformations (see Appendix |^, the 
prescription given in ( 2.12 ) indicates that the gauge invariant combinations are given hy Va '■= OeV^ (Ya — 4-4) 
and 11"^ := Ily/{a\^/io). The canonical first order Lagrangian for this degree of freedom is: 


= U^Va + K ( c|n 


tA/ 


1 , 4(™A + 2 + 2c|-H') 

(P)A + 2)"<"" - 


((2)A + 2) 


-2hEU^V, 


V^((2)A + 2)yA 

4 k CEa%yFJ 


(3.10) 


Solving for the velocity Va, we find the second order reduced action for 11^. It is convenient to express the diver¬ 
genceless vector as 11^ =: With this we obtain 


{v)g{ 2 ) ^ 2 k 


4Z(4) 


n 


= 2k 


(2) A 
(2)A-h2 
(2) A 


□ ■ 


2H' 

-,2 
E 


a 


HYx 


I (ai^n) I - 1 - ^) } («Bn) drjEdY, 


(3.11) 


where □ stands for the four dimensional d’Alembertian, is the d’Alembertian on the (2-1-1) dimensional de Sitter 
space, and K is the operator defined as 


K := -- 


9 + 


I 2 


dfjE'^ 2 

The above action is really very simple when expanded in eigenmodes. If one absorbs the factor 


2k(2)A 2K£{i+l) 1 

((2)A + 2) “ {e + 2){£-l) (a44’ 


(3.12) 


(3.13) 


by redefinition of 11, we basically obtain the action for an ordinary scalar field, n, living in the curved background 
describing the bubble geometry, with an rys-dependent mass term. 

We decompose the field n into modes as 

n = ^Aff7rp,„C/|l^'”(x), (3.14) 

pim 


where tv pim is the coefficient which represents the amplitude and t7|l^™(a;) is a suitably normalized mode function 
which is also an eigenfunction of the operators K and Then the renormalized field n is defined by IT := 

'^pim'^p^rnUE^ix)- i® Convenient to normalize U£^^{x) by means of the Klein-Gordon norm with respect to the 

renormalized field IT. With this choice of normalization, when we go to the quantum theory by setting the canonical 
commutation relation between the operator counter part of 11 and its conjugate, TVpim can be recognized as the 


anihilation operator which satisfies [TTpfm, = 

1. Writting in the form, 


q 

4 

II 

akvP{riE)y^^^{x% 

(3.15) 

the equation of motion separates into 



dS ^'yplm 

= (/ + 1)3^^^’", 

(3.16) 

k[vF] 

= vF. 

(3.17) 


These equations admit the following interpretation. The first tells us that behave as scalar fields of mass {p^ -1-1) 
living in a (2-1-1) dimensional de Sitter spacetime. The spectrum of masses is determined by (^.17), which is a one 
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dimensional Schodinger equation with effective potential Notice that the potential in (3.17) is positive 

definite and vanishes at infinity. Hence, the spectrum is purely continuous, with > 0. Therefore, our formalism 
bypasses the issue of a possible supercurvature mode discussed in ||h| . That mode was shown to be pure gauge, so it 
is not surprising that it does not arise with a suitable choice of variable which contains no gauge degree of freedom. 
If we choose the modes to be Klein-Gordon normalized in the (2+1) dimensional sense, i.e.. 


i J d^Xy^CE (^yp^-^ypi’^' - ypi’^ypi'm'^ = Su'dr, 


then the (3+1) normalization condition reduces to 

nOO 

/ u^uP' dr]E = S{p — p'), 


(3.18) 


(3.19) 


which is the standard one for the Schrodinger problem. For definiteness, here we choose 

ypem _ »^+^r(»^£+ 

where Vpe is defined by using the associated Legendre function of the first kind as 


(3.20) 


'Ppiip) ■= 


p: 


ip-4 


'■{ihE) 




(3.21) 


Notice th at th e factor (^^A +2 becomes zero when £ = 1. In this case, we have to go back to the original 
Lagrangian (HI)- From Eq. ( |3.8| ) we find H^ = 0, which means that one extra constraint arises. Therefore there 
remain no physical degrees of freedom for i = \ mode. In fact, this case can be quickly treated along the lines of 
Faddeev-Jackiw approach. Substituting A + 2 = 0 into the Lagrangian (3.7), and casting it into first order form, 
we obtain 


= u^Va + 


kce 


Ai 


U^UvA - 2hE VaU^ - Pe WaK 


(3.22) 


which is independent of Fa- The variation with respect to Wa gives the constraint Tly' = 0. Then the normalizability 
of the mode functions requires Hy = 0. After substituting this constraint, the equation of motion for Va becomes 
also a constraint, which enforces the Lagrangian for ^ = 1 to vanish. Modes with £ = Q are also absent from the action 
by construction. The absence of modes with £ = 0,1 is what we expect, because the odd parity mode represents one 
of the tensor degrees of freedom inside the lightcone, for which the modes £ = 0,1 do not exist. 

Now we relate the quantities in the outside of the lightcone with those in the inside of it. Inside the lightcone we 
can use the tensor harmonics to decompose the tensor part of the metric perturbation into modes. Thus the mode 
function defined in Appendix B |Q will be the most convenient choice of the variable to specify the tensor 

perturbation there. In order to relate the amplitude T^ptm to we compare the (pA)-component of the metric 

perturbation in the synchronous gauge (Vp = 0). Following the notation in Appendix D, we associate a subscript (or 
superscript) N to indicate the quantity evaluated in this gauge. From Eq. (3.8), is evaluated as 


h^'pA = -2 kceTIa 


= E 


2KCEAf-T^pem 

aE 




A||b 


u 


p 


(3.23) 


On the other hand, the explicit expression for (pA)-component of the tensor harmonics is given in |^,|T^. After the 
analytic continuation to region C, for the odd parity graviatational wave perturbation, we obtain 



{£ - l){£ + 2)T{ip + £ + l)r(-zp + £ + 1) 

2p2(p2 _|_ i)r(ip)r(—ip) 


CEVp,eA^Y^\^uljl 


(3.24) 


Hence we find that the amplitude TVpim 


is related to the variable U, 


(-) 


ptm 


by 
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1 


r{-ip + e + i) 


dU 


(-) 


y + i)r(ip +1 + i)r(zp)r(-zp) dijE 

Conversely, by using the equation satisfied by of ref |ll|] , 

1 d 


plm 


a| drjE 


a% 




u(-) is expressed in terms of tt as 


jA-) _ Ar:„-^+i /r(zp + £ +i)r(ip)r(-ip) _2_ d{aEuP) 


(p 2 + i)r(-zp + £ + 1 ) 




driE 


(3.25) 


(3.26) 


(3.27) 


B. Even parity modes 


The even parity modes contain two dynamical deg rees of freedom. One of them is the usual scalar mode, and the 
other is the even parity tensor mode discussed in [pl| . After lengthy algebra, complicated by the fact that the spatial 
sections are not homogeneous outside the lightconenom the nucleation event, the Lagrangian can be cast into second 
order form as the sum of a Lagrangian for the scalar mode plus a Lagrangian for the even parity tensor mode. The 
details of the reduction of the action are given in Appendix C. Here we only discuss the meaning of the final results. 

For the scalar part, we have 


5 C) ^ 1 j ,/Z:^(dq){^^n + :i-d]qdrjEd^x, 


where we have introduced the Schrodinger-like operator 




(3.28) 


(3.29) 


where stands as before for the d’Alembertian on the (2+1) dimensional de Sitter space of unit radius. The 
variable q is related to the gauge invariant potential $h of Bardeen Q when evolved to the outside of the lightcone 
(see Appendix BE 


q = -^^H. (3.30) 

K<Po 

Putting q = qeY^ dA'P q {x) with the mode function of the form = aYq^{VE)yP^"^{x^), the equation 

of motion separates into (|3.16|) and 


0[qP] = {p^+4)qP. 


(3.31) 


Just as in subsection III A, th e ma sses (p^ + 1) of the (2+1) dimensional fields are determined as the eigenvalues 

O =p^ +A=: (3.32) 


of the Schrodinger equation ( ^.31| ). Now, if we absorb the factor 

1 


(AAp )2 ’ 


by defining q := Y^pimUg we obtain the action for an ordinary scalar field. As before, we require that Ug ™ is 
normalized with respect to the Klein-Gordon norm for the renormalized field q. This normalization condition reduces 
to 


^^The potential is given by := —ip + {B — E') (see Appendix]^. We also recall that is eqnal to ^ of Kodama 
and Sasaki 0 and to — Qf Mukhanov, Feldmann and Brandenberger |^. 
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qPqP'dr]E = 


6(p-p') 

6ppf 


for continuous p, 
for discrete p. 


(3.33) 


Sinc e the potential of the operator O is not positive definite, we cannot determine the spectrum of p^ unless we solve 
Eq. ( 3.31 ). If the spectrum obtained previously by ignoring degrees of freedom of the metric perturbations does 
not change (except for the wall fluctuation mode at = —4; see below), the spectrum will be continuous for > 0 
and there may be one discrete mode at — 1 < p^ < 0. We will discuss this issue in a forthcoming paper [Q. 

Note that for the modes with = —4 we have 0\q\ = 0. If one ignores the metric perturbations, these correspond 
to the wall fluctuation modes H, la-M. Once the metric perturbations are taken into account, however, the wall 
fluctuation modes are found to be contained in the continuous spectrum of the gravitational wave perturbations 
Ip^ . Hence we expect the modes with p^ = —4 cease to contribute to physical fluctuations. In fact, there is a 
strong evidence that this is true by examining the regularity of the metric perturbations due to these modes [plf . 
Unfortunately, however, we have no rigorous proof for it. One possible (and probably most reasonable) stand point is 
to require the square integrability of the mode functions, otherwise integration by parts cannot be performed. Then 
we find that the discrete modes at p^ = —4 should be excluded from the spectrum. This can be seen as follows. 

For a while, we neglect a positive definite term in the operator O. Then Eq. (^.3l|) for p^ = —4 becomes 


OoM = 0 , 


dril 


^0 ■= “TUT +<^0 ( — 


The two independent solutions are easily found as 


rVE 


Ni N2 I , 2 , 

9i = — , 52 = — / Po d-qE, 

Po Po J -00 


(3.34) 


(3.35) 


where Ni and N 2 are constants. One readily sees that qi is badly divergent at qE = ±00 since Pq vanishes there 
(this is necessary for regularity of the instanton). As for 52 , it is regular at = —00 by construction but diverges 
at qE = 00 . Thus there is no square-integrable solution in the limit of weak gravitational coupling. Now, since the 
solution <72 has no node, the lowest eigenvalue pp + 4 of the operator Oq will be greater than zero. Then since the 
term we have neglected from O is positive definite, the lowest eigenvalue of Eq. ( 3.31 ) will be greater than p§ + 4, 


which is positive definite. Hence we conclude that no square-integrable solution of Eq. (3.31) exists at p^ = —4. 


If the operator O in front of q were absent from the action ( 3.28 ), the above fact would be sufficient to exclude the 
modes with p^ = —4. However, the action for these modes seems to be an indeterminacy of the form zero times infinity. 
On one hand, O vanishes, but on the other, q are not integrable on the spacelike surface x = 0. As we have mentioned, 
however, there are evidences that they do not contribute to physical fluctuations. Hence it seems reasonable to accept 
the square integrability as the guiding principle and exclude these special modes from the spectrum. 

For the tensor part, we have 


= 2k 




,dS 


(K w) 


dS 


a-1- K 
(2)A((2)A-k2) 


w dqE d^x 


(3.36) 


where K is the operator defined in Eq. ( 3.12 ). If we expand w by means of the eigenfunction of AT, the operator K 
can be replaced with the corresponding eigenvalue. Then we can absorb the factor 


2k K 


2k 


(2)A((2)A-k2) {i-i)e{e + i){i + 2) ■(Arf)2’ 


(3.37) 


by redefinition of variable and we obtain the action for an ordinary scalar field. As before w is decomposed as 
w = qe WpimU^^"^■ As in the case of odd parity, the spectrum is purely continuous, with p^ > 0. 

Now we relate Wpi^ with the mode function [|^ defined in the inside of the lightcone (See Appendix B). 

In order to relate the amplitude Wpim to 17^^, we focus on the traceless part of the (Ai?)-component of the metric 
perturbation, v, in the synchronous gauge {A = S = T = 0). As before, we associate a subscript (or superscript) N 
to indicate th e qua nt ity evalua ted in this gauge. 

From Eqs. (C17), (|D8[) and (D9) and with the aid of the equation 


( 2 )a(( 2 )A-h 2 ) 

2kKce\/w 


n„ 


(3.38) 
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which follows from the reduced Lagrangian (Cl^), ujy, is evaluated as 


KJyJ / ^ \ 

v'n = - aEcl(^)A((^)A + 2) ( 3 . 39 ) 

On the other hand, using the expression for the tensor harmonics given in @0, after the analytic continuation to 
region C, we have 


VN 


= -E' 


E 


pim 


12r{ip + i + i)T{-ip + i +1) 

+ i)r(ip)r(-zp) 


X {2hEdp + (-£(£ + 1) + 2hl - 2p^cl)) VpiY^^ul+l. 
Hence, we find that the variable w is related to the variable through the relation 

1 I 




r{-ip + £ + i) 


du 


(+) 


Y {p^ + l)T{ip A £A l)r(ip)r(-ip) dijE 


p£ 


(3.40) 


(3.41) 


As in the odd parity case, using the equation for which is the same as for \ Eq. (3.26), the inverse relation 
is given by 


U, 


(+) 


pi 




T{ip + £ + l)T{ip)T{-ip) _2 

(p 2 + l)r(-ip + £+l) 


d{aEuP) 

drjE 


(3.42) 


IV. SPECTRUM OF GRAVITY WAVES 


As an application of our formalism, let us find the spectrum of long wavelength tensor mod es pr edicted in open 
inflationary models. This reduces to solving the scattering problem for the Schrodinger equation ( |3.17 ). The potential 
for this problem vanishes at both rjE —> —oo and tje oo, so the asymptotic behavior at ±oo of the two orthogonal 
solutions for the energy p^ can be taken just as incident plane waves from ±oo with momentum p which interact 
with the potential and produce reflected waves which return to ±oo with reflection amplitude a±, and transmitted 
waves moving to Too with transmission amplitude g±. That is, in the limit r]E —> ±oo, the are given by 


P 

“(+) “ 




{g+e^P^^+e-^P^^) {^e ^+oo), 


cr+ e 


— iprjE 


iVE -oo). 


(4.1) 


and 


P 

M(_) - 




J_ gWE +Oo), 

{g_ + e^P^’^) {r]E ^ -oo). 


(4.2) 


Here we note that p is non-negative. Using the Wronskian relations we obtain 


|a+p = l-|p+|Y |a_|2 = l-|e_|Y 

CT_ =0-+, -I- aZp+ = 0. 


(4.3) 

(4.4) 


Using Eq. (4.4), we can show that 


drjE u 


P i/P - n 

(+)“(-) “ 


(4.5) 


and hence and are orthogonal. Note that normalization condition (3.19) is satisfied because 
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le±l 


1 ) = 1. 


(4.6) 


Analytically continuing the solution inside of the lightcone by means of tje = — 77 /? —i 7 r/ 2 , where rja is the conformal 
time in region R, the amplitude of perturbations well after the modes have crossed the horizon [rjn 0) is given by 


^+)l 




= — (coshTrp + 3? . 

TT 


(4.7) 


As we can see from the equation above, the bubble manifests itself in the spectrum through the real part of the 
reflection amplitude of the Schrodinger problem. 

In the thin-wall approximation, we can take the interior of the bubble as a pure de Sitter space, with scale facto r 
given by aE = coshry^;) outside the lightcone from N. In this limit, we can integrate out the potential in (3.17) 


and express it as a delta function with strength As = K/ii?vv/2, where /i = / drjE and Rw = ciE{ijw) are the 

surface tension and the radius of the wall, r espec tively. The reflection amplitude for this delta function potential is 
As/(2p + * As). Using Eq. (|3.42| ) with the fact that the scale factor inside the lightcone is given by 

a = — l/(i7sinhry/i), we recover the spectrum for given in reference [p^ : 




(+) 2 \ 


2TrKH‘^ 


^ p(p2-(-1) sinhTrp 


2kH^ coth TTp 

p(p2 + 1) 

where R is the reflection coefficient, given by 


1 - 


R 


cosh np 


■f-)t) 

As cos 2priu 


2p sin 2pp^ 


As 


R = 


(As)= 


4p2 + (As)^ 


(4.8) 


(4.9) 


Of course, the validity of Eq.(4.7) is not restricted to the thin wall regime, and a more complete analysis in the 
general case will be presented elsewhere Notice that this equation has been derived neglecting the term K,p'^(2 
in the equation for the evolution of the mooes inside the lightcone. This is justified for long wavelength modes, which 
are frozen in soon after bubble nucleation. Eor modes that enter the horizon at times t» H~^, which corresponds 
to p >> I, the generation of perturbations occurs during the second stage of inflation and has little to do with the 
bubble profile. Instead, the details of the slow roll potential will be important. 


V. CONCLUSIONS AND DISCUSSION 


In this paper we have applied Fadeev and Jackiw’s formalism for constrained systems to the problem of cosmological 
perturbations in the one bubble open inflationary universe (The cases of flat and closed universes have been also 


considered in Appendix 


)■ 


We have found the reduced action for a gauge invariant variable describing scalar degrees of freedom and for 
two gauge invariant variables describing tensor degrees of freedom. This tensor part coincides with the one found 
previously in ||ll|] . 

The nucleation of a bubble breaks the 0(4,1) symmetry of de Sitter space down to 0(3,1). It is known that, 
neglecting the self-gravity of the bubble, there is a special scalar mode with eigenvalue = —4 which corresponds 
to fluctuations of the bubble wall. This mode can be seen as the Goldstone mode associated with the breaking of 
symmetry 0(4,1) down to 0(3,1). We have seen that the wall fluctuation mode disappears from the spectrum of 
scalar perturbations once gravity is included. This is somewhat reminiscent of what happens in gauge theories: the 
Goldstone mode disappears when the gauge fields (gravitational degrees of freedom) are included (the tensor modes 
acquire a mass term on the wall, where the Goldstone used to live, which cuts off the infrared divergence encountered 
in j^). As pointed out in m (see also |3^), the disappearance of the wall fluctuation mode is perhaps not too 
surprising. Even in the absence of self-gravity, this mode can be written as a tensor mode, which contributes to 
microwave background anisotropies just like any gravitational wave would do The study of tensor modes in 

[ 0 ,H showed that in the weak gravity limit, the ‘infrared’ contribution of gravity waves to microwave anisotropies 
reproduces the effect of bubble wall fluctuations. 

As an application of our formalism, we have derived the spectrum of long wavelength tensor modes in open inflation. 
This is given in terms of the real part of the reflection coefficient in a one dimensional scattering problem, where the 
potential barrier is a function of the bubble profile. In the thin wall regime, we recover the results of |1^. A more 
complete study of this spectrum and that of scalar perturbations will be presented elsewhere [Q . 
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APPENDIX A: SECOND ORDER LAGRANGIAN 

In this appendix we derive the action for the fluctuations of a scalar field coupled to gravity on a FRW background, 
up to second order in the perturbation variables. 

The metric is written in the ADM form 

ds^ = —(A/"^ — NiN'^)drf' + 2J\fidx^dr] + ^^^gijdx^dx\ (Al) 

where M is the lapse, Mi is the shift function and gij is the metric on the constant rj space-like hypersurfaces in 
which we have foliated spacetime. Up to total derivative terms, the purely gravitational part of the action can be 
written as 

V = ^ =j^J - K\K^j) (A2) 

+ \d,{^/^gM^^^g^^)d, In ^^^g + d^Md^iV^g 

where 

d^ij = (A3) 

With a prime we denote a derivative with respect to time 77 , and the vertical bar stands for the covariant derivative 
with respect to the spatial metric ^^^gij- The action for the scalar matter field is 

Sm= J d*x^/^ • (A4) 

Now we expand the metric and the scalar field over an FRW-like background solution. The perturbed metric and 
the perturbed scalar field read 

ds^ = a{g)‘^{—{l + 2A)dg^ + 2Sidx^df] + ( 7 ^ -I- hij)dx^dx^, , 

‘P = Poiv) +dp, 

where 7 ^ is the metric on the constant curvature space sections. The background fields a and (po satisfy the equations 

n^-H' + JC = ^Po'^ 

2n'+H^+IC='^{-p'o^ + 2a^V{po)), (A 6 ) 

tpo" + 2 Hpo' + MV^{ipo) = 0, 

where Ti. := a'/a, and K, is the curvature parameter, which has the values 1 , 0 , -1 for closed, flat and open universes 
respectively. 

Expanding the total action, keeping terms of second order in perturbations, and using the background equations, 
we find 

S — Sgr + Sm =50-1- <525, (A7) 

where So is the action for the background solution and <525 is quadratic in perturbations: 

<525 = ^/ d^xa^^{-2{2n^ + 'H')A^ + ( 571 ,) - - (5*p - 

+ 2^{dip'^ — <5<^l*<5(/3|i — — 2k((^q< 5<^'A -I- a^U<^<5(/?A) 

+ 4(5*|, - - HA) + IC{2A^ - 2Ah\ - (A 8 ) 
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We have raised and lowered spatial indic es w ith 7 ^. 


In the open inflation case, comparing (A5) with (3.3), we see that the action for small perturbations outside the 


lightcone can be found just replacing a? 
the text. 


—a\ and 7 ,^ ^ expression (|A 8 |) . The result is equation (^) in 


APPENDIX B: REDUCTION INSIDE THE LIGHTCONE 

Although in the case of open inflation the t = const, surfaces are not Cauchy surfaces for the whole spacetime, it is 
known | 0 ] that they can be used to normalize the subcurvature modes, i.e. those modes for which the eigenvalue of 
the Laplacian on the hyperboloids of homogeneity and isotropy is smaller than -1. Furthermore, the resulting reduced 
action can be, in some heuristic sense to be discussed later, analytically continued to the outside of the lightcone. Then 
it is found that we obtain the correct result even for supercurvature modes, i.e. modes other than the subcurvature 
modes. Compared with the reduction in the outside of the lightcone, the analysis in the inside of the lightcone is very 
simple. Therefore, for an alternative less rigorous but rapid derivation, in this appendix we consider the reduction of 
the Lagrangian directly inside the lightcone. We shall simultaneously consider also the case of flat and closed spatial 
sections. —I —1 — 

Inside the lightcone the metric perturbations can be decomposed into scalar, vector and tensor modes il8| , |l9| , |2l| , 
regarding the way in which the modes transform under spatial coordinate transformations. On a homogeneous back¬ 
ground, the modes are decoupled in the action, and evolve independently. Thus we expand the metric perturbations 
as 




(Bl) 


where V'j B and E are scalar modes, Fi and Vi are vector modes,^ and tij is a tensor mode. Fj and Vi are divergenceless, 
and tij is transverse traceless (TT), i.e. 


F^li = V^l, = 0 = t\=f^y. 


Substituting the decomposition (Bl) in (A 8 ), the action is decoupled into three pieces 


(B2) 


62 S =(") <525 ^25 S 2 S. 


(B3) 


1. Scalar Perturbations 


The action for scalar perturbations reads 




— j d'^x - UnA^' + 2A^{2A - - 2{H' + 271^)A^ 

+k((5(^'^ + SipASip — a^VipipSip'^) + 2K{SipQ'tjj'Sip — A — a^V^^ASp) 
+/C(- 6 V^^ + 2A^ + Ut/jA + 2(B - F')A(B - F'))} 

+4A(B - F')(^(pgS<p -Ip' - HA), 


(B4) 


where A is the laplacian associated with 7 ^^. 

To apply FJ formalism we first have to cast the Lagrangian in first order form, defining momenta for the 
variables whose time derivative appears quadratically in the Lagrangian. As usual, the conjugate momenta are 
defined as 


Hw, := -4(")^25= ^ i-Zi}' + AE' + ?,'^ip'.5p>-AB-?,nA , 

S'lp K \ 2 / 

He := -^(^), 525 = (iCE' + 'ip' -- ICB+ nA) . 

SE K \ 2 / 


(B5) 


**It should be noted that Fi and Vi defined here are different from Fa and Va defined in Eq. (3.6). 
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The first order Lagrangian for scalar perturbations turns out to be 


(s)£ ^ _^(s) ^{s) _{s) 

= n^V'' + + ^eE', 

= 4»^/7(1 + 3K) (-*: nj + 2 n.n, + M + 2 (a + 3 

+ ((A+ 3/C)V'^ - I (^A + 3/C-H2-7f' + ^^ V) , 

Cb = n_E, 

Ca = H-(-(A + 3/C)^/> + - {H(Pq - ifio)Sip'j . 


^(p'o n^S(p 


(B6) 


We observe that neither A nor B enters into £i, so there is no dynamical evolution for these fields. They correspond 
to SA/"^ in the notation of the introduction. These fields appear linearly in the lagrangian, and their equations of 
motion, Ca/b = 0 contain no time derivatives. They allow us to evaluate two of the momenta in terms of the other 
fields. Moreover, the constraints Ca and Cb are the generators of the infinitesimal gauge transformations associated 
with diffeomorphisms. Under a scalar diffeomorphism generated by the vector = (A°,AI*), the metric transforms 
into g^i, + from which we can read the variation of all the scalar components of the metric perturbation. By 

commutation with Ca + XCb we recover the transformation law for the canonical fields: 

Sg'ijj = -nx°, (JgBv, = ^ (A + 3/C)A°, 

SgSip = (PoX°, SgHs^ = - (figB) X°, 

SgE = A, = 0. 


The constraints and the scalar Hamiltonian satisfy the following algebra: 

{('^^n,CA} = dgCA-ncA + CB, {^'^^kCb} = o, {Ca,Cb} = o. 


(B8) 


The time derivative of the constraints appears due to its explicit time dependence. This derivative acts only on 
background quantities but not on canonical coordinates. Under a gauge transformation, transforms as 


5g(")£i = -A° 


_d _ 

dg dg 


Ca-X 


_d _ 

dg dg 


Cb ■ 


(B9) 


Using these results and the fact that the action is invariant under a gauge transformation, we can recover the 
transformation law for the lagrange multipliers : 

5gA = X°' + nX°, SgB = X'-X°. (BIO) 


Now we proceed with the phase space reduction. We start by solving the constraints Cb = 0 and Ca = 0 for n_E and 
respectively. After substitution of the constraint, there is no E dependence in the Lagrangian, so the Lagrangian 

becomes a functional of only H^, if: and 6ip, ij), 5ip\. The dissaperance of E is related to the fact that 

there is no variable other than E itself whose gauge transformation depends on A besides the Lagrange multipliers. 
Hence it is not possible to construct a gauge invariant combination which contains E. Therefore E necessarily vanishes 
from (®)£*. .—^ 

Applying the formula given in (2.12) or equivalently looking at the gauge transformation law given in Eqs (|B7j), 


the gauge invariant combinations are found to be constructed from the remaining variables as 


T , ^ . 

= '0 H- 

■Co , 

2aV7 3. 

n^, = H^- ^(A- 


3IC)S(p. 


(Bll) 


The action expressed as a functional of II^ and can be obtained, up to total derivative terms, by simply putting 
(5(/3 = 0 in and replacing with 11^ and ■0 with ’S'. We finally obtain 
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- ^((A + sKjg- + ^Ti,f 




’S'(A + 3/C)’®' 


1 


lll^r 


k/C 


4a2^ "^(A + 3/C) 


n^. 


(B12) 


Notice that the procedure is equivalent to fixing a gauge where E and 5lp (the variables conjugate to the momenta 
we have solved the constraints for) are set to zero. 

For flat universes we can easily recover the results of or mi- If /C = 0, using the equation of motion for 
we can eliminate the momenta in favor of the velocity ’S'', and we will find the following second order Lagrangian: 


1 


= 


2 7^2 


which after the rescaling 




(’S''2 + ’i'A’4'), 


atpQ 


n ’ 


becomes the Lagrangian of a scalar field in fiat spacetime with time dependent mass, 


,-17’* 


/’( 2 ) _ 

- 2 


The reduced gauge invariant variable i? coincides with the one found in ref [p3[ : 


(B13) 


(B14) 


(B15) 


17 = a{Sip+^^). 

In the general case we can perform the canonical transformation 

2k77 1 


’S' = 


= 


_ 

4 aP-yf^Lp^ A + 3 /C 
2 




aV7'Pc 


0, 


r.n, - (A + 3 /C) q H -p, 
2 tL 


and solve for the momenta p to obtain the second order Lagrangian 


(2) _ 


— _ 
- 


(A + 3/C)q(n - m'^[a,ipo]) q, 


1 


i2[a, H ( 4/C - 2n' + <^0 


1 


(B16) 


(B17) 


(B18) 


Here □ stands for the four dimensional scalar d’Alembertian. The action for this lagrangian can also be written as 


J + 3^)gin + A + 3/C^ Qin dqdP X , 


^ / 2 / 


(B19) 


where <7in = a q, which is formally analogous to ( |3.28 ). Note, however, that (Oq) is changed by (A — 3)qi„. These 
operators are not the analytic continuation of one another, but both have the same eigenvalues on solutions of the 
equations of motion. 

We note that the reduced gauge invariant variable q is proportional to of Bardeen , 


9 = —r 

Kifo 


’S' - 


Hk 


2a2^(A + 3/C) 

For fiat universes (/C = 0), q is proportional to Q of ref p^ , 

Q 


U^]=^{i^-H{B- E')) =:-^ #H. 

' K(Pq Kipo 


(B20) 


(B21) 
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and the action (B18) reduces to the one found in ref. [p2| if we replace q with Q. If we take an ansatz of the form 
q = a-^ qP{r])Y^{x^), where are the scalar harmonics on the homogenious spatial section, the equation of 

motion separates into 


AYP^^ = (-p2 
d[qP] = (p2_4/C)qP. 


(B22) 


2. Vector and Tensor Perturbations 


As we have said, neither tensor nor vector modes couple to the scalar perturbations. The vector part of the action 
carries, as we will see, no dynamics. We find for it 


1 


(-)S2S = - — 
4k 




(B23) 


where 


Vm = V-A - 2/C Vm, Fm = ^-A - 2/C Fm. (B24) 

We can compute the corresponding first order Lagrangian, 

(B25) 

a Vt 

The field lAt has no conjugate momenta, and enters as a Lagrange multiplier which enforces Tf”* to vanish. After 
substituting tt™ = 0, we will end with a vanishing Lagrangian. In this model vector modes are pure gauge. 

For the tensor modes we find 

(‘^^25 =j-Jd^x a"V7 } . (B26) 

The action can be cast easily in a first order form, 

= J d\ ^ - 2/C)t,,)| . (B27) 

The Lagrangian is already in a canonical form and has no constraints. We can decompose Uj by using the normalized 
transverse-traceless tensor harmonics [^9|Jll| , and as 

= E ^ (B28) 

p£m p£m 

Then, it becomes manifest that there are two decomposed degrees of freedom for each p, i, to, which correspond to 
gravitational waves. 


APPENDIX C: EVEN PARITY PHASE SPACE REDUCTION OUTSIDE THE LIGHTCONE 


Inserting the mode decomposition ( p.q ) in the action (3.4), the even parity modes decouple from the odd parity 
ones. The even parity part reads 




+ 2C,') — 2 kp'q + ah. a — dli. w - — 


+ {w’ 

\ 2 , (-(2)A-2(2)Ahs2) r2 . 

S^-2S'w+^ -5-AJ-(2)a^'T 

J CE^ 


(2) A 

2 c £;2 


Ce 
(2)AT2 

2 c £;2 


+ T (2) A A -f 2 K (2) A A (2 + ( 2 ) a) v' + 2 A He C' 
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+ 2(2)AC'+ I -2/i|C" + 2c£;2 {-2nA + K‘p'^5^-w') c -%- 

+ C ((2 +(2) a) (w + (2)Aw)-2A (^(2)a + 2c£; 2^ -2(2)A/l£;^ + 4/l£;i 
-2 hE w) + (2) A e" - ^ (_2 (2) A HeA- (2) a (2 + (2) a) ^ 

+ 2(^)A A - + 4 CE^ (^nA - ^(fig S(fJ w' + ^ -hw + 4ce‘^'^ A 

A‘2, He A 2 v4 A — ((».A A 4 ce^'^ a ‘2 He 

+ 2hE^'^''AAv -\ - ^^“^2 —^ A 2c£;^ (l A 2H^ A H') A'^ — 4KaE^ S^p Ace^ 

— 2kSpce‘^ Ao ^ — a^ ce"^ Kvv') ~ A K(5a^| • 


By a overdot we denote a derivative with respect to p, and by a prime we denote a derivative with respect to rjE- 
The conjugate momenta are introduced as 


_ a^CE\AIj ^ feg(w A A 2'HS A A 2hECj 

a%CE^ (,0 - 2i + 2S' + 4HS + 2 (seC + A|) , 

<jkAAA + 2)i + 2hEA - (»)& + 2k) , 


II 7 ’ =- 


n,„ = -- 


n„ = 

= a^CE^/^ {Sp - p'oS) ■ 
The corresponding first order Lagrangian reads: 

(A£ = (A£^ _(A fi_CgS- Q5 - CqC, 

^ ^>Ci = 11^424 A IIt’T' a n^A A nt,v A n^^A, 
a\cE\/^ ( ^^^AT' 


where 


m = -- 


2k [ 


CE 


A T (-4 (2) A H 24 - ^ _^2k^‘^^A6p Aq 

\ qe^ ce Vw 


— ^^^A w' — A A 2 ) v'^ +w A A 4c_e^^ A A 


K He (n^ A 2 n^) 




(2)Au) 


[Hew) 


4 K /l£; n„ A 


0^2 Ce y/A 
A 4 c£;^ {h a — — Ao 


2^/2 


2 (2 A (^)A) a£;2 cb 4 /aJ 

a(2)Au ( f(2)AA4c£") Aa + 2n^) \ _ ^"^A(("^Aa2) ce 

\\ ) aE‘^R\Ao J 2 

/ie2(2)^%2 ^2242 ((2)AAc£;2 (2 A 2 (2a(2^A) 7i:2 a2 7i:'A^^^AH')) 


2 A (2)A 


2 /2 

_^C£_w- \- k6p‘^ - qe^ ce^ “ 2k6pce‘^ Ao “ 4«;a_E^ ^A Ace^ V, 


-kce' 


Ap’^ } 


a%CE\/io \ ( 2 )A 


CE^ Be" 


— Bu, — 


n,. 


Ba" b 


6tp 


(2)A(2 a(2)A) 


2 K 


Cs = AoBii^ — 2Hn„, — fly — qe 
Cj = B„ — A A CgBy, 


f A j' + 4A£iA A(„, + (»)&„) + (»)&„' + iArj 

Vue/ k V 4 / 


(Cl) 


(C2) 


17 



Cc = -2hE ( i((2) A + 2)(«; + (2) Az;) + 2cIHw' + clw" + 2^4 A' - (2) AA 


\2 


+ 2c|( 27^2 +7^'_ i)A- (2)AT' - Ka%c\V^^^5if - nc^ip'^d^p' ) . 


Notice that when £ = 0,1, equation (|C2| ) is meaningless because of the factors involving the laplacian. For the 
moment, we assume that £ ^ 0,1, and postpone the discussion of the fate of this two modes until the end of this 
section. 

We note that the helds C, and S do not appear in the canonical form Ci, and only appear linearly in the lagrangian, 
i.e., they are i5A/’^-like variables. The constraints are Cs = 0, Q = 0 and C(^ = 0, which generate even parity gauge 
transformations (see Appendix ^). 

To reduce the phase space we proceed following the way discussed in section |I| By taking linear combinations 
of the constraints Cg = 0, Q = 0 and Q = 0 we can construct which takes the for m Cg = — 'P^i[q\, where 

n^,IIiu} and {g} = {A,(5<p,T,u, w,n^,II t},. Now we can apply the formula ( 2.12 ) to obtain the gauge 
invariant combinations of variables as 


where 


$ 

T 


11$ 


Fix’ 


1 


Sp 
a-E —r 
^0 


A - 

ttE 

5p 

T+4-- 

V’o 

kHe 


2 / 
CeV , 


'hCc^w -\- c^w T 


(2) A 


-w 


fix 


^ 5p ^ w' 


kHe 


V4 Pq 


w = w + 2Ti.-^ + 
‘fio 


(C3) 


(C4) 


Substituting the constraints into the original first order Lagrangian to eliminate p^ = !!„, 11^}, we obtain the 

Lagrangian that depends only on q. Then, simply taking 5p = v = w = 0 and replacing A, T, n^i and fix by T, 
n$ and IIx, respectively, we finally get an action which depends only on the two pair of canonically conjugate gauge 
invariant fields: 


(s)£* = n$ # + Hx T n*, 

^ ^*sl _ 2hxn; ^ _ niif _ n$ 

2ax^cx-\/w (2)A((2)A + 2) (2) a (^(2) A + 2) ax^ cx-^0; ue'^ ce 

kce 2 (-At hE nx - A ax" ncEs/^^)T ax"4JT(^)AT 

(2)Aax^v^ 2kce 

ax^ ce ((")A + cx"(2 + ((^)A + 2){2U^ + H'))) # 

«:((2)A + 2) dax^cxVw' 

where If*, and 11^ are functions of T, IIx, 4* and 11$ determined from the constraints and given by 


= 

n: = 

n* = 

T = 


nx + n'^-Hn$ a|;/ixVw(2)AT 
4 k.cep'q 

„2 JT, I alcx^^^AyAl 

a%CE\fio 
2tihE 

(2) A$ + 2 c|(l - 21-e 2Uc 


a\'HcE\JTo ^ 
nhEPo 




I $' + (2)A(27i:T + T'). 


(C5) 


(C6) 


To disentangle the two degrees of freedom, we proceed in the following way. This Hamiltonian carries the counter¬ 
parts of the scalar and even parity tensor modes inside the lightcone. Since they are decoupled there, we can expect 
that they are also decoupled outside the lightcone. Hence we choose new coordinates {s, i;} such that when they are 
analytically continued inside the lightcone, they reduce to pure scalar and pure tensor variables in a particular gauge. 
A convenient choice is the longitudinal gauge for the scalar modes and the synchronous gauge for the tensor modes, 
in which Si = 0 and the traceless part of the spatial metric perturbation inside the lightcone becomes a purely tensor 
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quantity. Let us call this gauge the Newton gauge for convenience. Then it is easy to see that $ is a purely scalar 
type quantity. For tensor modes, we use the fact that Si — 0 hence T = 0 and u is a purely tensor type variable in 
the Newton gauge. Then with the help of the equations of motion, we can show that the following pair of fields s, v 
become purely scalar and purely tensor type variables (see Appendix ^ for a brief discussion about this subject): 


s 

V 


aEf'o 

K[aET] - . 

driE 


(C7) 


Then we can find new canonically conjugate momenta IIs, such that the Hamiltonian decouples into two 

pieces. Just for reference, we recall the definition of the operators O and AT, 


d = - 


dvl 





+ 2^0 


/ 2 


(C8) 


It is useful to keep in mind the following relation between O and K: 


^-d = k-^ 


(C9) 

dr]E <fo ipo dr]E 

To find the appropriate momenta basis in which the Hamiltonian decouples, we propose an ansatz for it, namely 


IIs + + Q[$] + 

Ut = a£;iF[n„]+f[T]+A^[$], 


(CIO) 


where Q and T are in principle arbitrary differential operators, but A# and Xt must be related in order to keep the 
t ran sf orma tion canonical. The momenta dependence of the transformation is found by requiring the transformation 
( |C^ - (|C10| ) to be canonical. Now we compute the canonical equations of motion for s and v using the old basis, and 
express the result in terms of the new basis. We find that s is independent of n.„, and that v is independent of Us. 
Then we define the operators involved in the definition of the new momenta basis in order to completely decouple 
these two equations. It can be shown that if we choose the operators as 


Q[$] = 
Xt[T] = 


A - 2(7f^ - l)c|) ^ 
kKe ’ 

a_ECgVcJ((^)A + 2c|) d ^ ^ ^ 

kHe drjE kHe 


(CII) 


s and V turn out to be 


aECEy/oj{^'^^^ + 2 c^) 


T[T] = — ~T, 




kHe 


kHece 


XM = 


cbv^C^^^A + 2c|;) d 2a%'Hc\^/u 


kHe 


drjE 




kHe 




1 


CeV^ 

K 

2 CE\ko 




4iFc|((2)A + 2-(l + iF)c|)\ p 

-(2)A((2)A + 2)- ) 

c|((2)^((2) A + 2) + 2(4 + (2)A(I - k))c\ - 4(1 + k)c%) 
(2)A((2)A + 2)/i£; 


K[v\ 


(C12) 


(CIS) 


which are already decoupled. It can be verified that the IIs and n„ defined by ( |C10| ) with the help of ( |Cll| )- (|C12| ) 
are canonical conjugates of s and v, so the two equations we have computed are two canonical equation of motion of 
the system. Therefore, to find the Hamiltonian in the new basis, we only need to know the two remaining canonical 
equations of motion. Computing IIs and 11^ we obtain: 
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o[tis 

K[U, 


— c|;V^ 


(2)A - ' 

(3 + — 5 —) s — 0[s] 


2 c%^K ^ (2)A(4+(^)A (1 - K)) - 4 c|(£( 2 )^ _ 2 + (1 + K)cl) 


(2)a((2)A + 2)/i| 


(")A + 34 

4k 


V — 


K^cl 


(C14) 


(2)A((2)A + 2)/li5 
_^(2 )a((2)a + 2) - 4(1 + i?)4) 


(2(4 + (2)A(l-i?))4 


Note that all the dependence has been absorbed in the differential operators O and K. Expanding ( C13| ) and 
(C14) in terms of eigenfunctions of this operators, we can read directly from them the coefficients of the Hamiltonian 
in the new basis. The corresponding first order lagrangian is 


*4* = 4[s,n3 


■ , TEi;] , 


(CIS) 


o 


n. 


Cq = IIsS - n* 

ZCEWOJ 


^ kK 

Trldi, — 

2 


'-E 


tJ 1 


o 


3 + 


(2) A 


-O s, 


1 - 


4lc|((2)A + 2-(l + if)c|)' 
(2)A((2)A + 2) 


n„ 


-U,,K 


c|((2)a((2) A + 2) + 2(4 + (2)A(1 - K))cl - 4(1 + K)c%) 


(2)a((2)A + 2)/ib 

1 ___ 2 c|Vw ( 2)^(4 (2) A(l _ K)) - 4c|((2) a£ - 2 + (1 + K)c%) 


2 K 


(2)A((2)A + 2)4 


(^)A + 3c| 
Kc% ^ 


The lagrangian Cg can be put easily in second order form. Solving for the momenta IIs, and defining q through 


s = 0[q], 


(C16) 


we find equation ( 3.28 ) of the text. 

The lagrangian needs a little extra work. Performing the following canonical transformation 


He 


n. 


(2)A + 2cI 


n. 


we find for 


ceV^^ kK 

KK{2{Kcl-hl)-(^)A) 
(2)A((2)A + 2) 


^ 2h|((^) A +c%)-K ((^) A + 2 c|)c| 
n, + 2c£;4J (2)A((2)A + 2)hB 


(C17) 


= n,„ w + Kwk 


C£4J((2)A-(l+i?)c|)^,. „ (2)^((2)^^2) 

W 11^, 


(2)A((2)A + 2) 


4:KKcE^/ijj 


-n. 


(CI 8 ) 


Solving for the momenta we find equation ( 3.3q ) in the texE _ 

When ^ = 0,1, simply looking at the definition of momenta (|Cl|), we can see that new constraints arises, due to 
the fact that some of them vanish. For ^ = 0, Ht and n„ become zero. In fact, the second order lagrangian 
for £ = 0 is independent of T and v. In this case we are left with a Lagrangian that depends only on three fields 
plus two lagrangian multipliers, therefore only contains one degree of freedom. Applying the Faddev-Jackiw 

formalism, the action for this degree of freedom turns out to be the one for the scalar degree of freedom, Sq , with 
£ = 0. Similarly, if £ = 1 the lagrangian is independent of n, so we have four fields and three lagrangian 

multipliers. As before, *-*^£^=1 only contains one degree of freedom. As expected, in this case we recover the action 
54 for £ = I. This is consistent with the fact that w represents one of the tensor degrees of freedom inside the 
lightcone, so it must be absent for £ = 0,1. On the other hand, q represents the scalar degree of freedom, so it must 
exists for all £. 
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APPENDIX D: GUESSING THE NEW VARIABLES 


As we have said, the guess for the variables which disentangle the lagrangian (|C5|) is motivated by their expression 
when evolved inside the lightcone in a particular gauge (Newton gauge, defined in Appendix]^, where we know they 
are purely scalar-type or purely tensor-type variables. To show this, we need to derive some useful and well known 
relations between the scalar potentials^ $h and #a- By a subscript (or superscript) N we indicate that the quantity 
is evaluated in the Newton gauge. 

First we consider the variables inside the lightcone. In the Newton gauge Bn and En vanish, so the constraint 
Cb = 0 reduces to 


V’w - + HAn = 0. 


(Dl) 


Substituting it into the definition of 11^, Eq. (B5), we hnd that this momentum also vanishes. 


n^ = o. 


(D2) 


Expressing the equations of motion for 11^; and which follow from Eq. (B12), in terms of ipN and 5(pN, and 
eliminating between them, we find 


1 Kj 

-{aipN)' - -<p'q5(Pn = 0 . 


Comparing with (Dl), we recover the well known relation 

($A —)An = 4’n{= —^h)- 


(D3) 


(D4) 


The analytic continuation of these relations (|D3D and ( |D4| ) to the outside of the lightcone does not change their form. 

Now, retur ning to the outside of the lightcone, we justify our choice of variables. In the Newton gauge, evaluating 
$ defined in (C3) with the aid of (|D3D and (D4), we hnd 


$ = 


1 


-O 


aE^Po 


2aE I 
—rVN 
^¥>0 . 


(D5) 


so €>, and therefore s and q, are already a purely scalar-type variable. Recalling the dehnition of q given by Eqs. (C7) 
and ( |C16| ):^ = 0[q\/{aEP^'o)^ we hnd 

2aEipN 


q = 




(D6) 


By means of the fact 4 >h = —tl’N inside the lightcone, we hnd Eq. ( ^.3C| ) in the text. 

To hnd the tensor-type variable, we use the fact that in the Netwton gauge vn is a purely tensor-type variable. 
The strategy is to construct a combination of $ and T proportional to vn- First, recalling that Tn = 0, we evaluate 
T dehned in (C2) in the Newton gauge as 


T = 


KaEP>o 


ji^{aE'4’N)' + Cb v'n- 


(D7) 


Then if we dehne 


T := qeT - TjiaEtpN)'= aEC%VN, (D8) 

we hnd this quantity becomes a purely tensor-type variable. Acting with AT, we hnally hnd the desired variable 

V := K[f] = K[aET] - A 

dq (fi'o Kip'o 

= k[aET] - ^(aE^), (D9) 

dq 


where we have used the relation 


^^The potential ^a is given by ^a A -|- {a{B — E'))'/a, as defined by Bardeen |2o| . 


21 

















APPENDIX E: GAUGE TRANSFORMATIONS OUTSIDE THE LIGHTCONE 


We can verify that under an even parity gauge transformation generated by A° Cs + A^Q+AQ, i.e. a diffeomorphism 
^ -vvith A^ = (A°, A^, the canonical scalar fields transform according to 

5gA = A°' + HAO, 5gT = -A° - c|A', 

5gw = -2nX° - 2hE\P - AA, 5gv = A, (El) 

5g5ip = v?o A°. 


The scalar constraints and the scalar Hamiltonian satisfy the following algebra: 


{(*)7f,Cs} = 5pCs-HCc, 


- clC's - hCc + Q, 

{Ca^Cjs} = 0. 


(E2) 


Notice that the parti al d erivative with respect to p acting on the constraints only affects background quantities. 
Using Eqs. (El) and (E2), and the fact that the action (C2) is invariant under a gauge transformation, we find the 
transformation law for the lagrangian multipliers as 


SgS = X° -clXP', Sg^ = X-XP, dgC = XP + hEXP + nX°. 


(E3) 


For the odd parity modes, under a diffeomorphism generated by A^ = (0,0, A"^), where A"^ is divergenceless; 
A^im = 0, the fields transform according to 


SgVA = -c|A' 


E'^A’ 


SqFa = -Aa ■ 


(E4) 


The algebra satisfied by the Hamiltonian and the constraint is: 




The action is invariant if we transform W as 


dgWA = Xa- 


(E5) 

(E6) 
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